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Memory effects are a key feature in the description of the dynamical systems governed by the 
generalized Langevin equation, which presents an exact reformulation of the equation of motion. 
A simple measure for the estimation of memory effects is introduced within the framework of this 
description. Numerical calculations of the suggested measure and the analysis of memory effects are 
also applied for various model physical systems as well as for the phenomena of "long time tails" 
and anomalous diffusion. 
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An abundant number of different problems involving 
molecular motions in condensed matters can be formulat- 
ed in terms of the generalized Langevin equation (GLE). 
The latter can be exactly received in the context of the 
Zwanzig-Mori projection operator technique pj, as well 
as within the framework of the recurrent relation ap- 
proach Q from the equation of motion. As known, this 
approach is used in the description of dynamical phenom- 
ena, ordinary and anomalous (such as anomalous diffu- 
sion) transport in physical, chemical and even biophysi- 
cal complex systems [l|, y] • 

One of the key features of the GLE is the fact that it 
contains an aftereffect function, termed a memory func- 
tion. If the memory function is delta-function correlated 
("a white noise"), the GLE is reduced to the ordinary 
Langevin equation corresponding to the system with- 
out memory, and the time correlation function (TCF) 
related to the momentum degree of freedom has sim- 
ple exponential relaxation. So, a formalism based on 
GLE (and/or likewise, on the corresponding generalized 
Fokker-Planck equation) inherently contains memory ef- 
fects, which characterize the system of interest 0] . 

Memory effects can appear in the velocity auto- 
correlation function (VACF) either through the presence 
of an oscillatory behavior or by means of slowly decreas- 
ing correlations. Thus, the memory character of the dy- 
namics of liquids related to influence of the cage formed 
by nearest surroundings on the particle movement (the 
so-called "cage effects") belongs to the first case, whereas 
dynamical processes related to anomalous diffusion and 
long time tails in the behavior of TCF's correspond to the 
second situation. It is necessary to note that in the sim- 
ilar treatment the notion of "memory" is not completely 
definite. In particular, it is found that non-exponential 
relaxation m can be observed in glasses, simple fluids 



and supercooled liquids 0, liquid crystals 0], plasmas 
H, frustrated lattice gases Q, proteins [HJ, disordered 
vortex lattice in superconductors [ijj- However, in con- 
cordance with GLE treatment, the process with memory 
effects constitutes any dynamical process, whose TCF de- 
cays, at long times, according to a non-exponential law. 
As a result, the following questions arise. To what degree 
is the studied process not Markovian? Correspondingly, 
how large is the difference between the relaxation with 
memory and an exponential relaxation? How strong are 
the space and time non-locality effects in the system? The 
situation would become more clear via the characteriza- 
tion by an adequate quantity ^3], which would provide 
a numerical measure of manifest memory effects. In this 
Letter we present such a measure: it is obtained by com- 
paring the time scales of the decay of the VACF and its 
corresponding memory function. 

Thus, the GLE of Mori-type can be written as 0,03]: 
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Mt(t-T)v a (T)dT + F(t). (1) 



where uj^ = {\Cv a \ 2 ) / (\v a \ 2 ) is the second frequency 
moment of VACF with the Liouville operator £, and 
the normalized memory function M\(t) can be related 
to the stochastic force F(t) by means of the fluctuation- 
dissipation theorem, (F(i)F(O)) = mk B Tcj^ 'Afi(t). 
Here fcs, T and m are the Boltzmann constant, the tem- 
perature and the mass of a particle, respectively. Multi- 
plying Eq. by w Q (0)/(|u Q (0)| 2 ) and performing an ap- 
propriate ensemble average (...), one obtains the mem- 
ory equation for the VACF, i.e.: 
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Ah{T)a(t - r)dr, (v a (0)F(t)) = 0, 
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here the VACF a(t) = (v a (0)v a (t)) / (v a (0) 2 } appears, 
which is related to the diffusion constant by the integral 
relation 0-0 = (k B T/m) J °° a{t)dt. 

Next we introduce a dimensionless parameter which 
characterizes memory effects for the relaxation process 
related to the velocity degree of freedom 
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where t 2 and are the squared characteristic relax- 
ation scales of the VACF and its corresponding memory 
function. Note that < 8 < oo. Obviously, the situation 
T a ^ t m corresponds to a memoryless behavior with a 
fast-decaying memory. Then, in accordance with Eq. (01 
we find that 8 — * oo. For processes exhibiting strong, pro- 
nounced memory effects, t 2 <§; r^, this parameter thus 
approaches 8 — > |l j |. 

Eqs. (J3J) can be rewritten in the terms of Laplace 
transforms of corresponding correlation functions, f(z) = 



Jo e zt f(t)dt, in following forms: 
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After a Laplace transform, Eq. Q reads: 
5(z) = [z + ^M^zf 1 , 
which allows one to obtain 

a'(z) + a(z) 2 



M[(z) 



uj^a(z) 2 



(5) 
(6) 

(7) 



where a'(z) = da(z)/dz and M{(z) = dMx(z)/dz. Then, 
taking into account the relations © we recast Eq. J^J as 



8 = co^ 



lim 



a'(z)a(zf 



o a'(z) + a(z) 5 
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Note that the knowledge of a(t) and its first derivative 
at long times, i.e. for z — > 0, and uj( 2 \ are sufficient to 
evaluate the measure 8. The complete knowledge of the 
memory function is thus not required. Moreover, the fre- 
quency moment of the VACF uj^ 2 \ and thus 5, is direct- 
ly related to physical observable quantities such as the 
distribution function and the potential of inter-particle 
interaction [lfij . 

Let us consider two contrasting situations: i) the sys- 
tem has a short-range memory, and ii) the system is char- 
acterized by physical (long-range) memory. 

In the first case (the memory-free limit) the memory 
function can be presented as M\(t) = 2rx8(t), Then Eq. 
<|2J is reduced to the ordinary Langevin equation (l^ : 



m^v a (t) + mLj {2) Tiv a {t) = F(t) 



(9) 



with a simple exponential solution for VACF: 

a(t) = e-" <2)T1 *, (10) 

which is correct for the VACF of a free Brownian particle 
with the relaxation time r = (uu^ti)^ 1 = m/7, m and 
7 are the mass and the friction coefficient, respectively. 
Let us define the parameter S for this case. Applying the 
Laplace transform to the memory function and VACF 
[Eq. tnH ] for this case, and taking into account Eqs. © 
one can find that t 2 — Tq and rj^ = 0. Then, Eq. (jHJ 
yields exactly 6 — > 00. 

The opposite situation is appropriate for systems that 
exhibit indefinitely large memory, i.e. the so-called strong 
memory limit. The time dependence for this model can 
be presented adequately as follows: 



Mi{t)=H(t) = 



1, t > 
0, t <0, 



(11) 



where H(t) is the step Heavisaid function. Introducing 
Eq. Illlll into the convolution integral of Eq. (0, one 
obtains the following equation 



da(t) 
dt 
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a.' 1 / a(r) dr, 



which has the simple solution 

a(t) = cos(\ 



(12) 



The similar VACF describes oscillations of the lattice 
particles in crystalline solids. Using the Laplace trans- 
forms of Eqs. 11H and 11211 we define from Eq. © the 
squares of characteristic relaxation scales of VACF and 
memory function, t 2 — l/uj^ and t 2 [ — > 00. Then, 
from definition (EJ we find 5 — ► 0. Thus, we have shown 
that the parameter 8 possesses true meaning for both the 
memory-free and the manifest, strong memory limit. 

We can therefore state with confidence that the quan- 
tity 8 sensitively detects and measures the presence of 
memory effects. The offered measure of memory can be 
obtained theoretically, experimentally and with the help 
of computer modeling if a(z) [or a(t)] is known. We next 
evaluate 8 for various model systems, thereby demon- 
strating its usefulness. 

A. The Rubin model The particle- impurity with 

the mass M is located in an infinite, ID harmonic lattice, 
particles of which have the same mass m (K is the cou- 
pling constant). The crucial quantity here is the mass 
ratio, q = M/m > 1. In case of q ^ 1 the movement 
of a particle with mass M is the same as the stochastic 
movement of a Brownian particle (the memory-free lim- 
it) P|. If the mass ratio q = 1, then the system consists 
of particles of the same mass, that corresponds to the 
model for acoustic phonons in crystalline solids at ordi- 
nary temperatures (the strong memory limit) Thus, 
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the Rubin model contains both memory-free and strong 
memory limits, transition to which is defined by value 
q. The term a{z) and the frequency parameter uj^ have 
the following forms here Q: 



a(z) 



(q - l)z + ^Jz 2 + AK/m 
From last expression we find 



(2) 



a(z = 0) 




a'(z = 0) 



M 



2K/M. (13) 



(q-1). (14) 



Then, from Eq. JSJ we find the expression 



*=|I ? -H 



(15) 



related 5 with the quantity q. It is seen from the given 
relation that for the strong memory limit at q = 1 we 
have 5 = 0, and for the memory-free one we obtain S — > 
oo. The squared g-dependence in Eq. ((15)) defines the 
transition in this model from one limit to the other. 

B. Ideal gas. TCF of density fluctuations for an ideal 
gas Il5j is 



(Ip(o)I 2 ) 



exp(-ft 2 t 2 /2). 



(16) 



Note that the Laplace transform of <fi(t) is related to 
the Laplace transform of autocorrelator for a longitudi- 
nal component of velocity of local density fluctuations, 

a/(t), by 



(z) = [ z + n 2 Mz)Y 



(17) 



The term a;(t) is similar to VACF and its frequency mo- 
ment a/ 2 ) is related to flf by a/ 2 ) = 20f [6|. Moreover, 
applying the Laplace transform to Eq. 11 fit and taking 
into account Eq. ifTTjl one finds 



ai(z) 



exp[-z 2 /(2ft 2 )] erfc" 



with the following low-frequency properties 
2 2- 



lim aAz) 

z— >0 



\im a'i{z) - 



Finally, one finds from Eq. ||Hl that 

s 4(77 - 2) 



7r(4 - 7r) 



1.69. 



(18) 



(19) 



It is the evident that the relaxation process related to 
fluctuations of a longitudinal velocity (momentum) com- 
ponent in an ideal gas is characterized by pronounced 
memory and independent from model parameters such 



as u> 



(2) 



C. Interacting 2D electron gas at long wavelengths and 
at T = 0. In Ref. Q it was shown that TCF of den- 
sity fluctuations and TCF of a longitudinal momentum 
component for this model are 



0(t) = J {2jn 2 t), ai{t) 



(20) 



where O J = k 2 e F , ep is the Fermi energy and k is mea- 
sured in units of the Fermi vector fop, Jn ls the Bessel 
function of the nth order. It is necessary to note that in 
the case of such time dependence of at(t), this TCF and 
its memory function have the same time behavior. As 
a result, they have the same relaxation time scales and 
6 = 1 that is easily obtained by the Laplace transform of 
Eq. <20ll with taking into account Eq. JjJJ, 

Long time tails. As known, TCF's in various physical 
systems can be characterized by long lasting tails jTflj l . 
For example, the long time behavior of VACF in liquids 



lim a(t) ~ t 2 . 

t — >oo 



(21) 



where d — 2 and 3 are two-dimensional and three- 
dimensional cases, correspondingly. In accordance with 
the Tauberian theorem, this is equivalent to 



lim a(z) ~ z 2 1 . 



(22) 



Then, one obtains the following low-frequency properties 



a{zf - 1, a'{z) = 0, 



if d = 2, 



a(z) 2 ~ z, a'(z) ~ —=. if d = 3, 
\ z 



(23) 



and from Eq. (jHJ we find that S = for both cases. It is 
evident, that the long time dynamics related to the frac- 
tal law ll2ll is characterized by memory effects which ap- 
pear directly by conserved long- living correlations. From 
this point of view, the source of such behavior is in the 
presence of arbitrary fragments of the regular motion in 
the system dynamics likely to be found in the billiard 
model [U]. 

Anomalous diffusion. Let us consider anomalous diffu- 
sion phenomena on an example of a free particle coupled 
to a fractal heat bath - the model used to study both 
sub- and superdiffusion. We don't consider the superdif- 
fusion region including ballistic limit with a > 2 [2lj | . 
The VACF for this model can be written as 



a(t) = £ 2 - Q (-7c^ ), 7a = 



ttAq 



mk bT sxn(an/ '2) 



(24) 



where E a (t) is the Mittag-Leffler function, Ao is the 
strength of the coupling |22|. Ordinary diffusion corre- 
sponds to a = 1, superdiffusion is associated with motion 
at 1 < a < 2, while subdiffusion occurs when < a < 1. 
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Eq. l]24t for a / 1 has an inverse power-law tail at long 
times 23] 



»(*) 



7a 



(25) 



One can readily deduce that for subdiffusion with a < 1 
the parameter S — > 0. For example, at a = 1/2 we obtain 
the same result as for a long time tail with d = 3 [see Eq. 
112311 ]. Let us next consider superdiffusion with 1 < a < 2. 
Then, for a = 3/2 we have 



lim a(zY 

z—>0 



lim a (z) 



(26) 



Taking into account Eq. 1(26)1 we find from Eq. (|HJ that 
S — > oo. As a result, although the VACF behavior in 
this case strongly deviates from the exponential relax- 
ation pattern, its memory function decays rather fast in 
comparison with the VACF, i.e. superdiffusion dynamics 
is similar to an ordinary Markovian process 

In conclusion, the notion of strong dynamical mem- 
ory is rather characteristic for the statistical properties 
of the time evolution of many systems, and as demon- 
strated with our study it depends on physical features 
of the underlying process. From the point of view of a 
recurrent relation method, properties of memory can be 
associated with the topology of a Hilbert space, where a 
dynamical variable is represented as a vector Q. Memory 
effects then occur both within a finite-dimensional space, 
and in an infinite-dimensional one, where the VACF is 
dissipative. On the basis of the GLE formalism we pro- 
pose a convenient measure of memory effects in dynam- 
ical systems. This quantity is based on a definition that 
allows one to extract in a tractable manner the role of 
memory behavior for the dynamics, and moreover, allows 
one to test the validity of memory-free approximations. 
The simplicity of the given measure makes its applica- 
tion possible in the analysis of various complex systems. 
Although the role of memory effects for chaotic process- 
es is not yet settled, the measure 6 may provide a useful 
guide in establishing an interrelation of manifest memory 
for such fundamental properties as chaotic behavior |25j 
and ergodicity [26] . 
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